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We study global entangling properties of the system of coupled kicked tops testing various hy- 
potheses and predictions concerning entanglement in quantum chaotic systems. In order to analyze 
the averaged initial entanglement production rate and the averaged asymptotic entanglement vari- 
ous ensembles of initial product states are evolved. Two different ensembles with natural probability 
distribution are considered: product states of independent spin-coherent states and product states of 
random states. It appears that the choice of either of these ensembles results in significantly different 
averaged entanglement behavior. We investigate also a relation between the averaged asymptotic 
entanglement and the mean entanglement of eigenvectors of the evolution operator. Lower bound 
on the averaged asymptotic entanglement is derived, expressed in terms of the eigenvector entan- 
glement. 
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I. INTRODUCTION 

Looking for quantum signatures of the classical tran- 
sition from regular to chaotic dynamics is the field of 
quantum chaos [l|, 0- Recently it was suggested [!, 0| 
that entanglement production in a quantum system can 
be a good indicator of the regular to chaotic transition 
in its classical counterpart. 

In both of the invoked studies [H, H| it was observed 
that the presence of chaos enhances the rate at which 
an initial product state is getting entangled. In the pa- 
per 0] Miller and Sarkar considered coupled quantum 
kicked tops. Single kicked top is a thoroughly studied 
model in the quantum chaos literature. Depending on the 
strength of kicks its classical dynamics is either regular or 
chaotic. Miller and Sarkar studied a system consisting of 
two identical kicked tops with an additional weak inter- 
action between them. The strength of kicks was chosen 
in such a way that the classical phase space of a single 
top was mixed (consisted both of regular tori and chaotic 
regions). Two tops were initially in a product state of 
two spin-coherent states. The reason for this choice is 
that spin-coherent states have a good classical limit, what 
gives a chance to relate the classical phase space picture 
to the quantum description. The spin-coherent state of 
the first top was chosen to lie in a chaotic region while the 
spin-coherent state of the second one was varied from a 
regular to a chaotic region. The system was then evolved 
and it was observed that the rate of entanglement in- 
crease (measured as the von Neumann entropy of the 
reduced density matrix) was higher when the second top 
was placed in a chaotic region. More quantitatively, it 
was shown that for different quantum initial states the 
rate of entanglement increase was proportional to the 
sum of two positive Lyapunov exponents calculated for 
the corresponding classical distribution of initial points. 
These results supported the claim about close relation 
between chaotic behavior of classical systems and entan- 



gling properties of their quantum versions. 

Further investigations revealed, however, that there is 
no such direct relation between chaos and entanglement 
[1 !, 0, i, H 01, EH- In particular it was observed that 
it is rather a specific time correlation function than the 
Lyapunov exponent itself that determines the entangle- 
ment production rate. Strong chaos is able to destroy 
the time correlations even on a very short time scale and 
thus it diminishes initial entanglement production rate. 

Apart from studying the entanglement production 
rate, which is a quantity calculated from a short-time 
behavior of the evolved state, one can also study asymp- 
totic properties of entanglement, i.e. those appearing in 
the long-time limit. The question whether the asymp- 
totic value of entanglement is related to chaos was posed 
in [!§]• If the coupling strength between kicked tops is 
large enough to observe the saturation of entanglement 
within given time, the asymptotic value of entanglement 
is higher the greater is the chaos parameter (the kick 
strength). This observation was made starting with an 
initial state in form of a product of two spin-coherent 
states placed in a region which for low kick strengths 
was regular, contained a part of a separatrix for stronger 
kicks, and eventually became chaotic for very strong 
kicks. 

In the present paper we also concentrate on the case 
of coupled kicked tops. In most of the previous stud- 
ies the attention was focused on analyzing how an initial 
product state of certain spin-coherent states evolves in 
time. Here we would like to address a more general prob- 
lem. Namely, how does the strength of chaos influence 
the global entangling properties of the evolution of the 
kicked tops. By global we mean properties not depend- 
ing on a specific choice of initial product state. In other 
words this is the question of whether the entangling ca- 
pabilities of the transformation depend on the strength 
of chaos. Global entangling properties were analyzed in 
the kicked Ising spin chain model Q, where evolution of 
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random states was investigated, and decrease of entan- 
glement production rate with the increase of chaos was 
observed. 

Another approach often used [13| in determining the 
entangling power of an operation consists of analyzing 
its eigenvectors. In the case of periodically driven sys- 
tems these are the eigenvectors of an unitary operator U 
corresponding to the one-period evolution. The degree 
of eigenvector entanglement is then regarded as an infor- 
mation about entangling properties of the evolution. We 
show that the information about entanglement of eigen- 
vectors does not, however, give the full picture of en- 
tanglement properties of the evolution, but only a rough 
estimate of the asymptotic behavior of entanglement of 
evolved states (see Sec lVI A[) . 

In order to discuss entangling properties of the evolu- 
tion of the kicked tops, we evolve not a single product 
state, but a whole ensemble of random product states, 
chosen uniformly with respect to action of SU(di) <E> 
SU(c?2) group, where d\, d 2 are the Hilbert space dimen- 
sions of the two tops. We calculate the averaged asymp- 
totic entanglement approached by evolved states, and the 
initial rate of entanglement production. All pure states 
are treated on equal footing here. Spin coherent states 
are as good as any other pure state. This way of aver- 
aging gives us information about the entangling power of 
evolution, as defined in [14j . 

Additionally we shall calculate the entanglement pro- 
duction when initial states are products of two spin- 
coherent states, where each spin-coherent state, param- 
eterized by two spherical angles 9, <f>, is taken indepen- 
dently from the ensemble with uniform probability den- 
sity on the unit sphere. The two ways of averaging give 
qualitatively different results. We shall analyze the re- 
sults with the help of the perturbative approach 0, |T(| 
and the analysis of the entanglement of eigenvectors. 

The observation that is worth mentioning here is that 
even in a very regular regime the entangling power is 
extremely high - higher than in chaotic cases - both when 
discussing the asymptotic entanglement behavior and the 
initial entanglement production rate. 

Our results are another step to reveal the relation be- 
tween chaotic vs. regular motion and entanglement pro- 
duction. 



II. ENTANGLEMENT 

Entanglement is a purely quantum phenomenon, di- 
viding states of a composite quantum systems into two 
classes: those which can be written as products of some 
states of the (two in this case) subsystems: 

I*) = |V>1> ® |V>2>, (2.1) 

called product states, and all other which cannot be writ- 
ten in the form (|2.1|) but instead involve a genuine, non- 



trivial linear combination 

l*> =X>J#«>® (2.2) 

ij 

The definition above applies only when one deals with 
pure states. The notion of entanglement for mixed states 
is more subtle [l5l ]. and it is very difficult in general to 
determine whether a given mixed state is entangled or 
not. In this paper the state of the composite quantum 
system we consider is always pure. 

For a given pure state \SSf) it is easy to check whether 
it is entangled or not. Observe that if the state is a 
product one, averaging over one of the subsystems of the 
corresponding density matrix p = gives a pure 

state of the remaining subsystem: 

Pi :=Tr 2 (p) = (2.3) 

whereas the same procedure applied to an entangled state 
produces necessarily a genuine mixed state. This obser- 
vation can be further quantified by calculating the linear 
entropy for the reduced density matrix 

S L :=l-Tr(pl), (2.4) 

which vanishes for product states and reaches the maxi- 
mal value: 

Sf** = \-\fd, (2.5) 

where d = mfn(<ii, d%) is the Hilbert space dimension 
of the smaller subsystem, for the " maximally entangled 
state" (by definition this is the state which reduces to 
the "maximally mixed" state of a subsystem - the latter 
is characterized by a diagonal density matrix with equal 
entries) . In this sense S is a measure of entanglement for 
pure states of a composite system. 

Another often used measure of entanglement for pure 
states is the von Neumann entropy of the reduced density 
matrix: 

5Vn = -Tr( Pl log 2Pl ). (2.6) 

For product states iSvN = 0, as the reduced density ma- 
trix pi is also pure, while for the maximally entangled 
states the von Neumann entropy of the reduced density 
matrix takes the highest value S"^ x = log 2 d. Contrary 
to the linear entropy, the above measure of entanglement 
has a nice operational meaning in terms of the number 
of maximally entangled states that can be distilled from 
a given number of non-maximally entangled states [l6| . 

Actually all quantities which do not increase under lo- 
cal operations (ie. operations acting separately in each 
subsystem) and classical communication quantify in some 
way the amount of entanglement present in a state. 
These in general are called entanglement monotones [l7| . 
Linear and von Neumann entropy discussed above are ex- 
amples of such entanglement monotones for pure states. 

In the following we shall use the linear entropy (|2.4p as 
the measure of entanglement. We choose this measure, 
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instead of the von Neumann entropy (|2 . 6[) . as linear en- 
tropy is easier to calculate and there is a perturbative for- 
mula for initial growth of linear entropy 0, E(| in weakly 
coupled systems which we shall use. Furthermore, in the 
investigations concerning relation between chaos and en- 
tanglement, where both von Neumann and linear entropy 
where calculated 0, EH, EH n0 qualitative difference in 
the behavior of the two was found, thus the choice of 
either of them is not crucial. 

It is argued that the presence of entanglement is impor- 
tant in many novel applications of quantum information 
processing [19J, which explains the prominence enjoyed 
by this phenomenon in many recent investigations. In 
our study we shall concentrate only on the interplay be- 
tween production of entanglement in a composite quan- 
tum system and its chaotic properties. 

III. COUPLED KICKED TOPS 

The kicked top is a paradigmatic model for stud yin g 
quantum chaos in finite-dimensional Hilbert spaces [201 ] . 
It is a particle with the total spin j and the dynamics 
generated by the Hamiltonian 

k °° 

H = pJ y + Y J 2 z ]T 6(t-n). (3.1) 

n— — oo 

Here J y and J z are the components of the angular mo- 
mentum operator fulfilling the standard commutation re- 
lations [J y , J z ] = iJ x , etc. The time dependence takes the 
form of an infinite train of delta-shaped pulses ("kicks") 
perturbing the free rotation periodically in time. The 
quantities p and k are adjustable parameters of the 
model. The latter, called the kick-strength, is scaled by 
the total spin j - observe that the total angular momen- 
tum J 2 = Jj + J 2 + J z is conserved, [J 2 , H] = 0, hence 
we can investigate the dynamics of the model for each 
value of j independently, restricting the discussion to 
the appropriate \{2j + 1) x (2j + l)]-dimensional space. 
To exhibit various interesting dynamical aspects of the 
model it is enough to change one of the parameters. In 
the following we put p = it/2 and vary k. 

The unitary time evolution operator transporting the 
wave function of the kicked top in time over one period 
of the perturbation, 

U = cxp (-i^i) ex P Jy) . (3-2) 

generates the Heisenberg equations of motion for the an- 
gular momentum operators J x , J y , and J z 

J' x = WJ X U =\{J z + iJ v ) e ^C^)(.4-i/2) + h. c . 5 
J' y = utj yU ^^(-j z + lJy ) e -^/ 3 )(J,-i/2) +hCj 

J' z = WJ X U=-J X , 

(3.3) 

giving the operators J x , J' yl and J' z at time t = n + 1 in 
terms of their predecessors J x , J y , and J z at time t = n. 



As in all studies of quantum chaotic phenomena we are 
ultimately interested in comparing quantum and classi- 
cal dynamics of the model. In the present the Planck 
constant has been set to unity hence the classical limit 
corresponds to j — ► oo limit (" large quantum numbers" ) . 
More formally one defines the quantities X, Y and Z as 
averages of J x / j, J y /j, and J z / j calculated in the initial 
state of the system. In comparing classical and quantum 
behavior it is reasonable to take as an initial state some 
minimum uncertainty state, in belief that (at least in the 
large j limit) the evolution of averages will be well repre- 
sented by the single classical trajectory starting from the 
point in the phase space around which the initial quan- 
tum state of minimal uncertainty is concentrated. Ob- 
viously, there is no particular reason to distinguish such 
states when investigating purely quantum properties like 
entanglement production. Appropriate minimum uncer- 
tainty states for spin j particles (so called angular mo- 
mentum coherent states) can be generated from the \j,j) 
state [i.e. the common eigenstate of J z and J 2 with the 
eigenvalues m = j and j(j + 1) respectively] by unitary 
rotations 

\9A) = (l + \l\Y jeliJx ~ iJy) \j>j)> 7 = e i0 tan^, 

(3.4) 

and the above described procedure leads to the following 
classical mapping: 

X 1 = Zcos(kX) + Ysva(kX), 

Y' = -Zsm(kX) + Ycos(kX), (3.5) 

Z' = -X. 

A detailed analysis of the classical dynamics (|3.5[) is given 
in [20(, let us only summarize that the system is inte- 
grable for k = and becomes visibly chaotic when k > 2. 
For k around 3 the phase space exhibits well developed 
mixed structure with few regular islands embedded in the 
chaotic see. When k w 6 islands of stability, although 
present, are very small and the chaos can be treated as 
fully developed for all practical purposes. From the point 
of view of quantum chaos, the islands are negligible if 
their phase-space area is smaller then the effective value 
of the Planck constant in our case), which will be 
the case in our calculations. 

In order to achieve our ultimate goal, ie. the investi- 
gation of parallels between chaos and entanglement we 
follow the idea of Miller and Sarkar Q and consider two 
coupled kicked tops with the Hamiltonian 

H = H 1 + H 2 + H I , (3.6) 

where Hi and Hi are the Hamiltonians of independent 
kicked tops (13. ip expressed in terms of the operators 
J Xl , J yi , J Zl and J X21 J V2 , J Z2 pertaining to each individ- 
ual top, whereas Hi is a nonlinear coupling term 

oo 

Hj = -J Zl J Z2 V 5{t-n). (3.7) 
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The procedure of obtaining the classical evolution equa- 
tions is exactly the same as the one for a single top de- 
scribed above, and yields [HI]: 

X[ = Z x cos A12 + Y\ sin A i2 
Y[ = -ZisinAia + Y1COSA12 

x\ = Z 2 cosA 21 +y 2 sinA 21 (3 ' 8) 
Y 2 ' = -Z 2 sinA 2 i +r 2 cosA 2 i 
Z 2 = -X-2, 

where Ai 2 = kX\ + eX 2 , A 2 i = kX 2 + eX±. In the most 
of the following the coupling strengths e will be small in 
comparison with k, it means that the degree of chaos in 
the system is determined solely by properties of dynamics 
of individual tops. 



initial product states. In order to treat all pure product 
states on equal footing one should choose the ensemble 
of product states with probability distribution invariant 
under the action of SU(di) x SU(d 2 ), where dx, d 2 are the 
Hilbert space dimensions for the subsystems. In this way 
we obtain a natural ensemble of random product states. 

It is not a surprise that SU(2) x SU(2) and SU(di) x 
SU(d 2 ) averages can lead to different quantitative esti- 
mates of the entangling power; what is more important 
they differ also qualitatively. 

One can also think about other characterizations of 
entangling capabilities of quantum evolution operators. 
For example, entanglement properties of eigenvectors of 
U can give some information about possible entangle- 
ment production. The matter, however, is rather subtle, 
as it will be clear from the subsequent discussion. 



IV. ENTANGLING POWER 

The main idea behind quantifying the entangling 
power of quantum evolution is to measure the ability to 
produce an entangled state out of an initial product state 
in the course of the quantum evolution. Although for 
particular reasons or applications we can choose a con- 
crete initial state and follow evolution of its entanglement 
properties when time passes, such a history would defi- 
nitely bear a lot of imprints of the initial state we chose 
to start with. In our study we are more interested in en- 
tanglement capabilities of the system itself, so it is more 
reasonable to take the average over some set of initial 
states - the idea advanced by Zanardi [l4[ 

e p {U) = {{S[U{\i )l )®\i> 2 ))))) M . (4.1) 

In the above formula S is some appropriate measure of 
entanglement (in our case it will be the linear entropy), 
and ((■}}tf >1 ,ip 2 denotes averaging over a set of initial prod- 
uct states IV'i) <8> \ip2)- 

The averaging procedure, however, needs more de- 
tailed reflection. As mentioned in the previous sec- 
tion, when investigating the quantum-classical corre- 
spondence, it is reasonable to take as an initial state 
a spin-coherent state and, consequently, a product of 
two such states for a composite system. The averag- 
ing amounts to integrating over the whole set of spin- 
coherent states parametrized by the spherical angles 9i 
and </>i, i = 1,2 in (|3.4[) . In order not to distinguish 
any particular initial state, we average over products of 
two spin-coherent states each parameterized by 6i and </>j 
taken independently from the ensemble with the uniform 
probability density on the unit sphere. This kind of av- 
eraging will be denoted as SU(2) x SU(2) averaging, as 
our ensemble is invariant under the action of rotation in 
either of subsystems. 

As already written, there is no particular reason for 
such a choice of the set of initial states when global en- 
tangling production properties of our system are inves- 
tigated. Instead one can average over the whole set of 



V. NUMERICAL RESULTS FOR COUPLED 
KICKED TOPS 

In this section we present the main results of the nu- 
merical calculations of the evolution of the coupled kicked 
tops. The spins of our two tops are respectively j\ = 19.5 
and j 2 = 20. This choice is an effect of a compro- 
mise. The spins should be high enough to allow clas- 
sical correspondence and low enough to be numerically 
tractable. Other authors, who considered higher spins 
(j = 40, j = 80), were able to perform their calcula- 
tions as they were evolving only a few different states. 
Our calculation are performed on ensembles consisting 
of several hundreds of states. Consequently we sacrifice 
the spin magnitude for the sake of being able to perform 
averaging over many states. The chosen spins ji, j 2 are 
not equal. This does not change evolution significantly 
(as compared with ji = j 2 = 20), but removes degener- 
acy among eigenstates of one-period evolution operator 
U. Lack of degeneracy is essential for the eigenvectors 
entanglement analysis, which will be explained later. 

The coupling strength is chosen to be e = 0.01, while 
the strength of kicks k (equal for both tops) will be varied 
from k = to k = 6. Chaos enters the classical dynamics 
of a single top at k w 2. The coupling constant is small 
enough to assure that all chaotic behavior is due to the 
kicks of the tops and not their interaction. 



A. Entanglement evolution 

In Fig. [T] and Fig. [2] we show calculations for 1000 iter- 
ations of tops evolution. In Fig. Q] the evolution of linear 
entropy averaged over 100 random initial spin-coherent 
product states is presented. 

During the first 1000 kicks most of the curves saturate 
to some asymptotic value. The only exception, k = 0.01, 
requires a little more time to saturate. The saturation 
of entanglement in the evolution of the kicked tops is 
present also in the case of single spin-coherent product 
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FIG. 1: The time evolution of the linear entropy averaged 
over initial spin-coherent product states (S'su(2)), calculated 
for different chaoticity parameters k of subsystems (differ- 
ent kick strengths). Spin magnitudes of two interacting tops 
are taken respectively ji = 19.5, ji = 20 and the coupling 
strength e = 0.01. The evolution time comprises 1000 kicks. 
The statistical limit (SL) for entanglement is drawn with the 
solid line and corresponds to the average entanglement of a 
randomly chosen state of the composite system. Increasing 
the chaoticity parameter causes in general an increase in the 
asymptotic value of entanglement (with the exception of the 
case k = 0.01). Initial entanglement growth is extremely slow 
for very regular dynamics (k = 0.01) 



state evolution - the averaging over spin-coherent states 
is not necessary, yet the averaged curves have smoother 
behavior. Different curves saturate to different asymp- 
totic values. This observation was made in [l|, where 
it was pointed out that higher is the chaos parameter k, 
the higher is also the asymptotic value of entanglement 
reached by the spin-coherent product state \9, </>) <g> \9, </>) 
(9 = 0.89, cf> = 0.63). Our averaged results confirm this 
observation only partially. While the asymptotic value 
indeed increases with k for k > 1, it is also quite high 
in the nonchaotic regime - k = 0.01. For k = 0.01 the 
asymptotic value is higher than that for k = 1, k = 2. 
Opposite to the nonmonotonic behavior of the asymp- 
totic values, the initial entanglement production rates 
seem to increase monotonically with k. We shall discuss 
these observation more thoroughly in the following. The 
statistical limit (SL) calculated as an average entangle- 
ment of a random pure state of the full system [21| is 
drawn with a thin solid line. 

In Fig. [J] we present the result of averaging over 100 
random product states with the probability distribution 
invariant under the action of the SU(di) x SU(g?2) group, 
where d\ = 1j\ + l,e?2 = 2j2 + 1 are the dimensions 
of subsystems (the state of one top is chosen indepen- 
dently from the state of the other). We generated ran- 
dom product states using random unitary matrices, dis- 
tributed uniformly according to the SU(di) x SU(d2) Haar 
measure, applied to a fixed product state (22J. Differ- 
ent asymptotic values of entanglement for different k is 
again visible. A monotonic increase of the asymptotic 



FIG. 2: The time evolution of the linear entropy averaged 
over the ensemble of initial random product states (Ssv(d)), 
with the probability distribution invariant under the action of 
the SU(di) x SU(d 2 ) group, where di = 2ji + 1, d 2 = 2j 2 + 1 
are the dimensions of the subsystems. Different curves corre- 
spond to different chaoticity parameters k of the subsystems 
(different kick strengths). Spin magnitudes of two interacting 
tops are taken respectively j\ = 19.5, ji = 20 and the coupling 
strength e = 0.01. The evolution time comprises 1000 kicks. 
The statistical limit (SL) for entanglement is drawn with the 
solid line and corresponds to the average entanglement of a 
randomly chosen state of the composite system. The initial 
entanglement growth is higher for regular dynamics (low k) 
than for the chaotic one (high k). While asymptotic values of 
entanglement are high for chaotic cases, the asymptotic value 
for the most regular dynamics k — 0.01 outperforms all other 
cases. 



values can be observed for k = 2,4, 6. However, the dif- 
ferences in the asymptotic values for k — 4, 6 are tiny. 
The k = 0.01 case is especially interesting. During 1000 
kicks the entanglement saturates to an extremely high 
value - higher than for all other values of k. Again, a 
nonmonotonic /c-dependence of the asymptotic entangle- 
ment is observed. 

Studying the initial production rate of entanglement 
one can observe anticorrelation with the parameter k. 
The fastest initial growth of entanglement corresponds to 
k = 0.01, and the slowest to highly chaotic cases k — 4, 6 
- chaos suppresses the initial entanglement production 
rate. 

For shortening the notation we shall denote the 
SU(2) x SU(2) averaging by the SU(2) averaging and the 
SU(di) x SU(d 2 ) by the SU(tf) averaging. The S\J(d) aver- 
aged behavior of initial entanglement growth is strikingly 
different as compared with the SU(2) averaging. In the 
latter case the initial entanglement production rate was 
almost zero for low values of k, while in the former it was 
extremely high. Consequently, one should always distin- 
guish between the entangling power of an evolution and 
its particular entangling properties in acting on a certain 
group of states as these two may behave very differently. 

Summing up the qualitative discussion we conclude 
that the entangling power of the evolution [correspond- 
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FIG. 3: This figure presents the dependence of asymptotic val- 
ues of entanglement on the chaoticity parameter k. Asymp- 
totic values of entanglement corresponding to averaging over 
initial random product states [the SU(d) averaging] are de- 
noted by stars and asymptotic values corresponding to the 
averaging over initial spin-coherent product states [the SU(2) 
averaging] are denoted by triangles. The dashed line repre- 
sents the Lyapunov exponent of the classical dynamics of a 
single top. The mean entanglement of eigenvectors of the 
evolution is denoted by rectangles. Statistical limit (SL) for 
entanglement is drawn with the solid line and corresponds to 
average entanglement of a randomly chosen state of the com- 
posite system. Chaos indeed increases asymptotic value of 
entanglement in the case of spin-coherent states and, to some 
extent, also in the case of random product states; neverthe- 
less very regular regime k ~ also manifests high asymptotic 
entanglement which can be related to high entanglement of 
eigenvectors. 



ing to SU(g?) averaging] both in terms of the asymptotic 
value and the initial growth is the highest for very low k 
- i.e. for very regular dynamics. 



B. Asymptotic behavior 

We give here more detailed results on asymptotic en- 
tanglement. In Fig. [3] the asymptotic values of entangle- 
ment are presented for different values of k. For more 
credible results the asymptotic values were obtained as 
averages of the linear entropy over the evolution of the 
tops between 50,000 and 100,000 kicks. Stars correspond 
to SU(d) averaging while triangles correspond to SU(2) 
averaging. The Lyapunov exponent obtained from the 
classical dynamics of a single top is shown by the dashed 
line. For the sake of later discussion we also included 
the mean entanglements of eigenvectors of the evolution 
operator - denoted by rectangles. 

In the case of SU(2) averaging, for very low values of 
k, the asymptotic values are high. With the increase of 
k they decrease, reaching minimum for k w 1.5, subse- 
quently, with the onset of chaos, they increase again and 
eventually saturate. The saturation value is a little below 
the statistical limit (SL). 

The SU(g?) averaging reveals almost no dependence of 



FIG. 4: Dependence of initial entanglement production rate 
on the chaoticity parameter k. Results obtained after averag- 
ing over initial random product states [the SU(d) averaging] 
are denoted by stars and results due to averaging over ini- 
tial spin-coherent product states [the SU(2) averaging] are 
denoted by triangles. The dashed line represents the rescaled 
Lyapunov exponent of the classical dynamics of a single top. 
In the case of random product states the highest initial en- 
tanglement growth corresponds to very regular dynamics and 
diminishes with the increase of the chaoticity parameter k 
while in the case of spin-coherent product states it is very 
regular dynamics that has the slowest initial entanglement 
growth. 



the asymptotic value of entanglement on k. Nevertheless, 
there is also a tiny dip around k w 2 and the entangle- 
ment for very low k is a little bit higher than that for the 
strongly chaotic regime. 

The behavior of the mean eigenvector entanglement, 
also reveals a minimum around k ss 2. Remarkably the 
entanglement of eigenvectors for very high k is signifi- 
cantly smaller than asymptotic value of entanglement in 
this regime. 



C. Initial behavior 

In order to grasp quantitatively the initial behavior 
of entanglement, we fit a line to the points representing 
short-time entanglement produced for certain value of k. 
Although the character of initial entanglement growth 
is linear only in the chaotic regime, while in the regu- 
lar regime it is quadratic (see Sec. IVI CI) , we perform 
the fitting in all regimes. This gives us sensible estimate 
of initial entanglement growth. The fitting is done for 
points corresponding to first 15 kicks. Regression coeffi- 
cients obtained in this way are shown in Fig. [4] for both 
SU(2) and SU(d) averaging, together with the rescaled 
Lyapunov exponent of a single top. For high values of 
k both averaging methods give the same results, which 
is caused by strong chaos which even during the first 15 
kicks is able to turn spin-coherent states into completely 
random states. For low values of k the entanglement 
production rate obtained with the SU(d) averaging out- 
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performs that of the SU(2) averaging. Furthermore it is 
approximately four times higher than in the strong chaos 
regime. We shall now move on to explaining these obser- 
vations. 



VI. THEORETICAL ANALYSIS OF 
NUMERICAL RESULTS 

Before moving on to separate discussions on the 
asymptotic and initial behavior of the entanglement pro- 
duction, let us first make a general remark which applies 
to both cases. Looking at Fig.[3]and Fig.[|]one can notice 
that when the chaoticity parameter is large (k > 5) dif- 
ferences between the SU(2) averaging and the SU(d) av- 
eraging disappear, both in results for the asymptotic and 
the initial behavior. It means that there is no difference, 
whether we choose as our initial states the ensemble of 
spin-coherent product states or the ensemble of random 
product states. Disappearance of the difference between 
two averaging methods for high k is due to strong chaos 
which very quickly turns initial spin-coherent states into 
random pure states. The distinction between two aver- 
aging methods, however, is crucial in mixed and regular 
regimes. 



A. Entanglement of eigenvectors vs asymptotic 
values of entanglement of evolved states 

Abstracting for a while from the case of the kicked 
tops, let U be a unitary operation acting in a Hilbcrt 
space Hi (g> H2, where dimensions of Hi, H2 are respec- 
tively d\ and d%. We denote by |ej), fa for i = 1 . . .did% 
the eigenvectors and eigenphases of the operator U. We 
assume here that the spectrum is nondegenerate. The re- 
duced density matrices of the eigenvectors after tracing 
out the second subsystem are pi — Tr 2 (|e,)(e.j|). As in- 
formation about entanglement of an eigenvector we use 
again the linear entropy of its reduced density matrix: 
Si = 1 — Tr(p|). The mean entanglement of eigenvectors 
is given by 

Seigen = 1-^-^1^). (6.1) 

i 

Assume that the initial state of the system is After 
n iterations of the operation U the resulting state reads 

\ip( n )) = y^exp(m<fo)(ej|y>)|e^). (6.2) 

i 

The reduced density matrix of the first subsystem corre- 
sponding to this state is given by 

p( n ) = ^exp[m(<^ ~ feW'K^le^ThfleiXejl), 
ij 

(6.3) 



the linear entropy of which reads 

S p (n) = 1 - Tr [p(n) 2 ] = 1 - Tr( E exp[in(fa - fa + fa - fa)} 

\jkt 

(e#)(^|e,)(e fc |^(^|e i )Tr 2 (|e l )(e,|)Tr 2 (|e fc )(e i |)). 

(6.4) 

For the asymptotic behavior of the linear entropy the 
important terms in the above formula are those which 
survive after the time-averaging. These are the terms in 
which phase factors disappear. Other terms oscillate and 
thus vanish when averaged over time. The nonvanishing 
terms correspond to either i — j,k = lari = l, k=j. 
Asymptotic value of linear entropy reads thus 

Sasymp - 1 - Tr ^£ | | ^) | 2 \ 2 PiP A - 

Tr (e ; |(e i |V>| 2 Ke J |V)| 2, &2(|e i )(e J |)Tr 2 (| ej )(e i |)^ . 

(6.5) 

We would like to stress here that the above formula is 
the result of the averaging over time of the linear entropy 
itself. One could take a different approach and first aver- 
age over time the density matrix and than calculate the 
linear entropy of the resulting density matrix. Following 
the second approach one would not get the third term 
in the above formula. This difference is important as it 
is certainly not the same to calculate the time-averaged 
entanglement of the evolved state or to calculate the en- 
tanglement of the time-averaged state. We argue that the 
first approach reveals better the entangling properties of 
the system evolution. One must admit, however, that in 
many cases the second approach will give qualitatively 
similar results. 

We would like now to relate the above formula for the 
asymptotic entanglement of an evolved state with the 
entanglement of eigenvectors. We start by proving the 
following inequalities: 

2Tr(Tr 2 (|e i )(e J -|)Tr 2 (|e i )(e i |)) < Tr(p 2 ) + Tr(p 2 X6.6a) 
2Tr(p iPj ) < Tr(p 2 ) + Tr(p 2 )(6.6b) 

Proof. We shall prove the first inequality. Let us de- 
compose the eigenvector |e,-) in a product basis |ej) = 
E c n in2 \ n i) ® K)- We can write: Tr 2 (|ei) (e 3 \) = 

C H C j1i . It follows that 

2Tr[Tr 2 (|e i )<e J -|)Tr 2 (|e j )<e i |)] = 
= 2Tr(C ra C J ' t C i C it ) = 2Tr(C nt C n C jt C J ) < 
< Ti'(C^C 2 C^C l ) + TriC^&C^C 1 ) = Tr(p 2 ) + Tr(p 2 ). 

In the derivation above we have made use of the Cauchy- 
Schwarz inequality. Assume that A and B are hermitian 
operators. The Cauchy-Schwarz inequality than reads 
Tr(^ 2 )Tr(B 2 ) > [Tt(AB)} 2 . Taking into account that 

[Tr(A 2 ) + Tr(£ 2 )] 2 > 4Tr(yl 2 )Tr(B 2 ), one arrives at: 
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Tr(A 2 ) + Tr(_B 2 ) > 2Tr (AB). This explains the deriva- 
tion above and gives (|6.6a[) . The proof of the inequality 
(|6.6bj) is analogous. 

We now return to the formula (|6 . 5[) . Using the inequal- 
ities (16.61) we can write: 



1 - Tr 



Tr ^|(e l |^)| 2 |(e J |^)| 2 Tr 2 (| ei )<e J |)Tr 2 (|e J )( e ,|) > 



>1-Tr [^{eiM^m 



PiPj 



mating entangling power from entanglement of eigen- 
vectors can be even more misleading. As an exam- 
ple consider a local transformation U = U\ (£> U2, 
which obviously has the zero entangling power. Let 
l/i) ® I/2) and \gi) ® | (72) be two degenerate product 
eigenstates. However, eigenstates in this situation can 
equally well be taken as (l/\/2) (|/i) <g> |/ 2 ) + \g%) ® \g 2 )) 
and (l/\/2) (|/i) ® I/2) — ® I52)), which are entan- 
gled states. Consequently calculating numerically eigen- 
states and their entanglement we could arrive at the 
conclusion that the entanglement of eigenvectors is high 
while the transformation itself has the zero entangling 
power. Nondegenerate spectrum is thus an indispensable 
condition in studying relations between entangling power 
of a transformation and entanglement of its eigenvectors. 



Tr K]|(e^)| 2 |<e J #)| 2 Tr 2 (| ei )<e J |)Tr 2 (|e J )( ei |) 



> 



> 



1 - 2]T |(e^)| 2 Tr(p 2 ) = 2^ |<e#>| a Si - 1, 



where Si is the linear entropy of the eigenstate | e^) . 

We are interested in averaging over initial states 
which are either spin-coherent product states or random 
product states. Both averaging procedures lead to 



((\{eiM 2 )),=l/d 1 d 2 



(6.7) 



Finally we arrive at the following relation between mean 
asymptotic value of entanglement of evolved states and 
mean entanglement of eigenvectors: 



^asymp _ 



1. 



(6.8) 



This inequality is valid for both the SU(2) averaging and 
the SU(d) averaging. Actually, as the equation (|6.7[) is 
valid for any ensemble of states which give a resolution 
of identity, so is the equation (|6.8p valid in all these cases 
and not only in the case of SU(2) or SU(d) averaging. 
This inequality puts the lower bound on the mean asymp- 
totic entanglement of evolved states, thus a high entan- 
glement of eigenvectors induce a high mean asymptotic 
entanglement of evolved states. 

One is not entitled, however, to claim that a low mean 
entanglement of eigenvectors implies a low mean asymp- 
totic entanglement of evolved states (compare Fig. [3]). 
To stress once again the inadequate information about 
entangling properties of a transformation obtained from 
studying its eigenvector entanglement, recall that even 
when all eigenvectors are product states, a transforma- 
tion can have a non-zero entangling power (for example 
the controlled-phase gate). 

This is a good place to comment on our decision to 
take ji different from ji- If we took identical spins 
for the two tops the degeneracy of eigenvectors would 
appear. The presence of degeneracy invalidate the in- 
equality (|6.8p . In the presence of degeneracy esti- 



B. Asymptotic behavior 

When explaining the asymptotic value of entanglement 
of an evolved state in a system of two strongly chaotic 
interacting subsystems one can formulate some statisti- 
cal predictions. One expects that states evolved in such 
systems tend to generically random states. This means 
that writing an asymptotic state, in a product basis: 



ni 



n2) 



(6.9) 



one expects the coefficients a ni „ 2 to be independent ran- 
dom variables. Depending on the symmetries of the sys- 
tem, one imposes different restrictions on the coefficients 
arriving thus at different state ensembles. It is then pos- 
sible to calculate the linear entropy averaged over such 
ensembles. [H[2l| 

Choosing an ensemble of states which are uniformly 
distributed with respect to the action of SU(<ii<i 2 ) group, 
one arrives at averaged linear entropy [21(: 



S, 



SL 



1 



d\d 2 



1 



(6.10) 



Consequently, the above expression is equal to the aver- 
age entanglement (measured as linear entropy of reduced 
density matrix) of a generically random state of the com- 
posite system. We shall use this expression, which we call 
statistical limit (SL), as a reference for the asymptotic 
entanglement of evolved states. 

When dimensions of subspaces are similar d\ ~ d 2 ~ d 
and large the statistical limit can be written in a simple 
form Ssl ~ 1 — 2/d. Compare this with the limit imposed 
solely by the dimensionality of the subspaces: S max = 
1 — 1 /d, where d is the dimension of the smaller subspace. 

In Fig. [5] we compare this formula with the values of 
the asymptotic entanglement for strong chaos k = 6, for 
different spin magnitudes of the two tops. There is no 
difference between values obtained for the SU(2) or the 
SU(d) averaging, and both are a little below the statisti- 
cal limit (SL). Notice that the eigenvectors entanglement 
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FIG. 5: Asymptotic entanglement of evolved states averaged 
over initial random product states [the SU(d) averaging, stars] 
and over initial spin-coherent product states [the SU(2) aver- 
aging, triangles) in the strong chaos regime k — 6, calculated 
for different spin magnitudes j. Spins of the tops are chosen 
ji = j and j2 = j + 1/2. Rectangles denote values of mean 
entanglement of eigenvectors. The statistical limit (SL) for 
entanglement is drawn with the solid line and corresponds 
to average entanglement of a randomly chosen state of the 
composite system. A big discrepancy is visible between the 
asymptotic entanglement of evolved states and the mean en- 
tanglement of eigenvectors. 
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FIG. 6: Asymptotic entanglement of evolved states in the case 
of averaging over initial random product states [SU(d) aver- 
aging, stars] and in the case of averaging over initial spin- 
coherent product states [SU(2) averaging, triangles] in the 
regular regime k = 0.01, calculated for different spin magni- 
tudes j. Spins of the tops are chosen: j\ = j and j'2 = j + 1/2. 
Rectangles denote values of mean entanglement of eigenvec- 
tors. Statistical limit (SL) for entanglement is drawn with 
a solid line and corresponds to average entanglement of a 
randomly chosen state of the composite system. Thanks to 
inequality (|6.8|l high entanglement of eigenvectors is an in- 
dicator of high asymptotic entanglement for both averaging 
procedures 



is hardly of any use here in predictions of the asymptotic 
entanglement. 

We noticed earlier that for low values of k asymptotic 
values of entanglement for the SU(2) and the SU(ci) av- 
eraging are also very high. In the light of the inequality 
(16. 8p this is due to the very high entanglement of the 
eigenvectors of the evolution operator. Asymptotic en- 
tanglements for k = 0.01 are presented in Fig. [^together 



with the entanglement of eigenvectors for different spin 
magnitudes. The statistical limit (SL) (|6.10[) is also plot- 
ted for comparison. 

The agreement between the statistical limit (SL) and 
the asymptotic values of entanglement of the SU(d) aver- 
aging is purely accidental. This underlines, however, that 
the regular regime can be as good (or better) in generat- 
ing high asymptotic values of entanglement for random 
states compared as the strongly chaotic regime. This 
agrees with the results in [23|, where different approach 
was taken - no specific model of dynamics was analyzed 
but instead the Hamiltonians of the subsystems where 
random matrices chosen either from chaotic or regular 
ensembles. 

It remains to be explained what is the reason for a 
very high entanglement of eigenvectors for low values of 
k. This effect of high entanglement of eigenvectors in 
the regular motion regime was also observed in the spin- 
kicked rotor system [241 ] and a mechanism is similar here. 
We shall give a qualitative explanation of this effect. For 
decoupled kicked tops e = 0, eigenvectors are product 
states. For low values of k there are many states that are 
nearly degenerate. Additionally, these product states are 
very close to the product states of eigenvectors of the J y 
operator. Indeed for k = one time step of the evolution 
of the tops is described by the unitary operator: 



Uq = cxp(— ipJ yi — ipJ t 



.1/2 I 



(6.11) 



so eigenvectors can be chosen as 1 2/1,2/2), the prod- 
uct states of the eigenvectors of J y operator. There 
are many degenerate eigenvectors of Uq, for example 
I - 3,j), I - 3 + l,i - 1), \j, -j) " these are 2j + 1 
states with the eigenphase (for simplicity we have as- 
sumed ji = j'2 and the value of the spin to be in- 
teger). Additionally there will be also groups consist- 
ing of 2j, 2j — 1, . . . , 1 degenerate states. The coupling 
term exp(— e J z \ J Z 2), expanded to the first order in e, 
will couple the state 12/1,2/2) with the states of the form 
\yi -1,2/2+1), |2/i + l,2/2-l), |2/i + l,2/2+l), 12/1-1,2/2-1) 
(as J z is a sum of lowering and rising operators in the 
basis of J y eigenvectors) . A weak coupling between non- 
degenerate states will not cause much change in the form 
of eigenvectors. Let \vi) be the set of degenerate eigenvec- 
tors of Uq. Due to the coupling, new eigenvectors will be 
obtained by diagonalizing a matrix of the approximately 
following form: 



p e 
e p e 
e p 



Eigenvectors of such a matrix have large contributions 
from many different vectors Vi, in our case it means that 
in every subspace of degenerate product eigenvectors new 
highly entangled eigenvectors will emerge due to the weak 
coupling. 
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C. Initial behavior 

The initial entanglement growth in the system of cou- 
pled kicked tops can be well understood with the help of 
the perturbation theory developed in @, [13] ■ 

It should be clarified here that, although the two pa- 
pers 0, fiot contains similar results, the motivation and 
the scope of their work is a bit different. The main mo- 
tivation of the work by Fujisaki et al. 0] is to study 
entanglement production in weakly coupled chaotic sys- 
tems and as a model they consider the coupled kicked 
tops. This is also the approach we take in this article. On 
the other hand, Znidaric and Prosen consider a problem 
of stability of quantum dynamics of a composite system 
with respect to weak interaction between two subsystems. 
They consider different quantities reflecting the stability 
of quantum motion such as fidelity, reduced fidelity and 
purity fidelity, and analyze their behavior under regular 
or chaotic dynamics for different times of evolution (com- 
prehensive study of the problem of stability of quantum 
dynamics can be found in see also 

& 0,123). Be- 
cause the linear entropy (5), which is used as a measure 
of entanglement in Q and in the present work, is related 
to the purity fidelity (Fp) by the formula 5=1 — Fp : 
there is a close relation between results obtained in [7[ 
and [3. 

The essential quantity for understanding initial entan- 
glement growth is the product of the time correlation 
functions of uncoupled subsystems: 



D(n, m) — Ciin, m)C2(n, m), 



(6.12) 



where the correlation function of an uncoupled subsystem 
is defined 

Ci(n,m) = (zi(n)zi(m)) - (zi(n)){zi(m), (6.13) 

and Zi(m) = J Zi (m)/ji. The Heisenberg picture is used 
here and J Zi (n) denotes the operator of projection of the 
angular momentum on the z axis at time n for the sub- 
system i (after the n-th kick of uncoupled evolution). 

The perturbation formula for the initial behavior of 
the linear entropy reads 



5 perturb (t) = 



t t 

EE 

m— 1 



-D(n, m). 



(6.14) 



For strong chaos one can neglect the terms with n =/= m, 
as any correlation is quickly washed out. For n = m 
due to the chaotic character of the motion D(n, n) w 1/9 
and does not depend on n (more detailed analysis of 
the behavior of the correlation function can be found in 
@, 0, H El). Finally for strong chaos one obtains that 
the initial growth of linear entropy is linear, and the rate 
of the growth is given by: 
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FIG. 7: Initial entanglement growth rates averaged over ini- 
tial random product states [the SU(d) averaging - stars] and 
initial entanglement growth rates averaged over initial spin- 
coherent product states [the SU(2) averaging - triangles]. The 
results are obtained in strongly chaotic regime k = 6 and 
for different spin magnitudes j. The spins of the tops differ 
by 1/2: j± = j and j2 — j + 1/2. The prediction of per- 
turbative formula (|6.15jl is drawn with the solid line. Due 
to strong chaos there is no significant difference between the 
spin-coherent and the random states. 



This formula is valid for any initial state, actually not 
from the very beginning of the evolution, but only after 
a short relaxation time when the linear growth of the 
linear entropy emerges. In our case the linear growth 
appears very quickly, after first few kicks. 

The comparison of the above formula with the initial 
entanglement growth rate corresponding to the SU(2) 
and the SU(rf) averaging, in the case of strong chaos k = 6 
is shown in Fig. [7] The agreement is good and, due to 
strong chaos, there is no significant difference between 
two ensembles. 

In the regular regime, the formula (|6.14|) . yields initial 
quadratic increase of linear entropy, due to nonvanishing 
correlations C(n,m) for n ^ m [7], Q. In this regime 
there is qualitatively different behavior of linear entropy 
increase between SU(2) and SU(e2) averaging (see Fig.f?]). 

In order to explain the difference we have to estimate 
values of the correlation functions C(n, m) averaged over 
the SU(2) and SU(eQ ensembles of initial states. Notice 
the following formulas for a particle with spin j: 



I 



dHsu(2)(ip\Jz\ip) 2 
d^sv(d){iAJl\^) 
d^sv(d) (i>\Jz\ip) 2 



jV + 1) 

3 : 

t 
3 ' 

JV + 1) 

3 : 

j 
6' 



(6.16a) 
(6.16b) 
(6.16c) 
(6.16d) 



dt 



' i ■ ■ 



(6.15) 



where d = 2j + 1 is the dimension of the particle's Hilbert 
space. Proofs of these formulas are given in the appendix. 
The same equations can be written when J z is replaced 
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by either J y or J x . If, however, (ip\J^\ip) is replaced by 
a mixed term - for example (tp\J z J x \ip) - the integrals 
(|6.16a| . (|6. 16c|) vanish. The same happens if in the inte- 
grals (|6.16b[) . (|6.16b[) the term (tp\J ' z \ip} 2 is replaced by 
a mixed term - for example {i})\J z \ip) {tp\J x \^>) . 

In the most regular case (k = 0.01) the evolution of 
a single top is mostly determined by its free rotation 
around y axis. Each period corresponds to a 7r/2 rota- 
tion. Consequently in the Heisenberg picture operator J z 
evolves approximately in the following way: J z (l) ~ J x , 
J z {2) « - J z , .7,(3) « - J x , .7,(4) » J z . . .. The value of 
the correlation at equal times C(n, n) is just the disper- 
sion of J z (n) operator. In the case of the SU(2) averag- 
ing this will be proportional to the difference of integrals 
(|6.16a[) and (|6.16b[) and thus proportional to In con- 
trast, C(n,n) for the SU(ci) averaging is proportional to 
the difference of integrals (|6.16c|) and (|6.16d[) and thus 
proportional to 1 - the same order of magnitude as in 
the chaotic case. The same fact will hold for correlation 
functions calculated for times differing by even number of 
kicks, whereas correlation functions calculated for times 
differing by odd number of kicks will be zero. As D(n, m) 
is a product of two correlation functions for the two tops, 
it will be approximately jxjz times greater for the SU(rf) 
averaging than for the SU(2) averaging. In the case we 
consider, this amounts to the ratio of initial entanglement 
production rate of approximately 400. This explains the 
strikingly different behavior of the initial entanglement 
growth for two considered kinds of averaging in for the 
k = 0.01 case. Different character of entanglement pro- 
duction rate for coherent and random states was also 
noticed in [9(. 

SU(g0 averaging reveals decrease in the initial entan- 
glement growth rate with the increase of chaoticity pa- 
rameter k As the states we average over the ensemble 
of random states invariant under the action of SU(rf), 
the increase of parameter k does not change the value of 
the dispersion of J z operator; thus C(n,n) does not de- 
pend on k. At the same time the increase of k decreases 
correlation for different times C(n, m) (n ^ m), and con- 
sequently the increase of the chaoticity parameter checks 
initial entanglement growth. 

In the case of SU(2) averaging, consequences of the 
increase of the chaoticity parameter are twofold. On 
the one hand the increase of chaos decreases correlations 
C(n, m) for n ^ m, and in this way transforms initial 
entanglement growth character from quadratic to linear. 
On the other hand it drives well-localized spin-coherent 
states into more smeared states with larger dispersion of 
J z operator, which results in the increase of C(n, n) func- 
tion calculated at equal times. Depending on the specific 
values of parameters of the system: ji, j2, e, one of theses 
factors may bare grater significance, and it may happen, 
that the increase of chaos corresponds to either the in- 
crease or decrease of initial entanglement growth (see for 
example numerical results in 12] for different values of e 
parameter) . We chose the parameters that reveal the in- 
crease of entanglement growth with growing k for SU(2) 



averaging in order to contrast this with the behavior of 
entanglement growth in the case of SU(d) averaging, as 
for the latter case the increase of chaos always diminishes 
initial entanglement growth. 



VII. CONCLUSIONS 

In this paper, the problem of the interplay between 
entanglement production in a quantum system and its 
chaotic properties was analyzed using the model of the 
coupled kicked tops. 

Entangling properties of the coupled kicked tops sys- 
tem were investigated by observing the evolution of 
two different ensembles of product states. Consider- 
ations of ensembles consisting cither of product spin- 
coherent states [SU(2) averaging] or random product 
states [SU(d] averaging) lead to qualitatively different re- 
sults in terms of the initial entanglement production rate 
and the asymptotic entanglement of evolved states. The 
asymptotic values of entanglement are high in a strongly 
chaotic regime (high kick strength), still they are even 
higher for a very regular regime (low kick strength) , and 
reach minimum for the parameters of the evolution that 
correspond to the onset of chaos. 

The SU(2) averaging reveals strong dependence of the 
asymptotic value of entanglement on the kick strength, 
while asymptotic values obtained after the SU(rf) aver- 
aging are quite insensitive to variation of this parame- 
ter. This is due to more classical nature of spin-coherent 
states, which "feel" more directly the transition form reg- 
ular to chaotic dynamics in the classical motion. In the 
case of the SU(d) averaging, the classical transition from 
regular to chaotic motion has only minor influence on the 
asymptotic value of entanglement. 

In order to gain a deeper understanding of the prob- 
lems, the asymptotic entanglement behavior was related 
to the mean entanglement of eigenvectors of the transfor- 
mation. We have proved an inequality relating these two 
quantities. We have also pointed out cases where this 
quantities differ strongly and one is not entitled to infer 
anything about one of them from the value of the other. 

Finally we have studied the averaged initial entangle- 
ment production rate. The striking difference between 
the SU(2) and the SU(e2) averaging was explained with 
the help of the perturbation theory 0, [13] • In the regu- 
lar regime spin-coherent product states become entangled 
much slower than random product states, due to a sig- 
nificantly smaller value of averaged variances of angular 
momentum components. 

Regular to chaotic transition can indeed be observed in 
the global entangling properties of the kicked tops. These 
manifestations, however, are different depending on the 
ensembles of states considered. When considering ran- 
dom product states the increase of the chaoticity param- 
eter always diminishes initial entanglement growth. For 
spin-coherent states increase of chaos results in two com- 
peting tendencies: chaos drives coherent states into more 
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delocalized states and thus helps entanglement growth, 
but on the other hand it destroys time correlations in sub- 
systems, which checks entanglement growth. We chose 
parameters of the evolution in which first tendency pre- 
vailed and thus contrasted the behavior of spin-coherent 
vs. random states. 

When discussing asymptotic entanglement we have ob- 
served that both purely regular and strongly chaotic 
regimes enjoy a very high asymptotic entanglement. Rea- 
sons for high asymptotic entanglement in regular and 
chaotic cases are different. In the chaotic case, chaotic 
dynamics in subsystems allows for a coupling of arbitrary 
states and consequently an initial product state can be- 
come highly entangled as the dynamics is able to pene- 
trate Hubert space of the composite system with almost 
no constraints. In the regular case evolution of a state 
is much more restricted, but in the system of coupled 
kicked tops, due to a very high entanglement of eigen- 
vectors (see inequality 16. 8 j) the asymptotic entanglement 
is also very high. 
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APPENDIX: DERIVATION OF FORMULAS (t67l6j) 



Formulas (|6.16a[ I6.16c[) can be explained using sym- 
metry arguments. If instead of J|, we took J 2 = 
Jj + Jy + Jf, the mean value (ip\ J \i>) = j(j + 1), for 
any state. Both SU(2) and SU(d) averaging do not dis- 
tinguish any direction in space. Consequently averages 
of J 2 , Jy and J 2 should be equal. In order to sum up to 
j(j + 1) for J 2 operator each of these averages must be 
equal j(j + l)/3 

Let us now prove formula (|6.16b[) . A spin-coherent 
state can be obtained as a rotation of \j) state: \6, <fi) — 
R(8,(f))\j). This allows us to write 



dpsv(2)(ip\Jz\ipf 



4tt 



sin(0)(0,</>|J z |M> 2 



a rn \m), where a m = x m + iy m are complex coeffi- 

m=—j 

cients. The SU(d) averaging over |?/>) results in a random 
distribution of a m coefficients with the probability distri- 
bution: 



P(a_j, . . . , cij) = P(x^j,y^, . . . , Xj,yj) 



-J y sm(0)sin(0)|-j) s 



1 

47T 



sin(60cos 2 (60j 2 = ^ 



T s(i - J2 4, - E 

m=—j rn——j 



Marginal distributions can be calculated from the above 
formula. In what follows we shall need the two lowest 
marginal distributions of P i.e. Pi{x) and P2(x,y) [27]]: 



P 2 ( X ,y) = Ql { l- X *-y2fj-\ 



Note also the following equalities: 



dxdy x y P 2 (x,y) 



8(j + 1)(2j + 1) 
1 

8(j + 1)(2j + 1) 



(A. la) 
(A.lb) 



We are now prepared to prove formula (|6.16cl[) : 



da_j . . . dajP(a_j, . . . ctj) 



= —j d6dct>sm(6){-j\R(6,ct>)U z R(0,ct>)\-j) 2 = 
= — j d9d<f>sm(9){-j\J z cos(6>) + J x sin(0) cos(^) 



E a n{ n \ J * E a ™l r ' 
n=— j m=-j 

r ( 3 \ 

/ da,-j . . .da,jP(a,-j, . . .dj) I |a m | 2 m 

\m=-j J 



/ '} 3 
da-j . . .dajP(a-j, . . .a,-) E E l am l 2 l a ™ 



m——j n——j 



The proof of formula (|6.16d[) requires a little more ef- 
fort. An arbitrary state can be written as: \ip) = only the terms with either n = m or n = —m contribute, 
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as all others cancel out: 



f da-j . . . dcijP(a-j , . . . cij) |a m | 4 m 

m=—j 

3 

E 



a m \ 2 \a- m \ 2 m 2 = 
m——j 

dx-jdy-j . . . dx j dy 1 P{x^. j ,y^ h Xj,yj) 
j 

E ( x ™ + y 4 m + 2x2 m y 2 m ) m2 - 

m=—j 

j 

E { X m + ym)( X -m + V-mW = 



2 J dxx A P 1 (x)-2 J dxdy x 2 y 2 P 2 (x,y)j ^ 



3 

2 



m 

i(i + l)(2i + l) _ j 



2(j + l)(2j + l) 3 6- 
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